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seems controversial. 
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1. Ordering in two dimensions 



In the context of phase transitions, two-dimensional models exhibit a very rich 
variety of typical behaviours, ranging from conventional temperature-driven second 
order phase transitions (e.g. Ising model) to first-order ones (e.g. q > 4— state Potts 
model), with also specific properties of models with continuous global symmetry 
which may present defect-mediated topological phase transitions (e.g. XY model) 
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or even no transition at all (e.g. Heisenberg model). Models of nematic- isotropic 
orientational phase transitions belong to this latter category of systems displaying 
a continuous symmetry. Ordering in low dimensional systems is likely to be frus- 
trated by the strentgh of fluctuations. On qualitative grounds, consider e.g. how 
fluctuations develop within the framework of Landau theory when the tempera- 
ture decreases from the high temperature paramagnetic phase toward the transi- 
tion temperature. The response to a localized magnetic field applied at the origin, 
h5(x), follows from Ginzburg-Landau functional minimization of the free energy 
F[m(x)} = J cbc(|am 2 (x) + |6m 4 (x) + K\Vm(x)\ 2 — M(x)m(x)). After linearization 
and Fourier transform, it yields 

r 2 m(k) + |k| 2 m(k) = A (i) 

where £~ 2 = a/2K is the correlation length. The response m(k) is here proportional 
to the correlation function Fourier transform G(k) and the fluctuations are measured 
through 

k B T X = V-\{M 2 ) - (M) 2 ) = £G(k). (2) 

k 

The latter sum is converted to an integral from to some cutoff A. It is diverging 
with f in Id (£ k G(k) ~ £Arctan (Af)) and in 2d (E k G(k) ~ hi(A£)) while it is 
bounded in 3d (^2kG(\a) ~ const.). The divergence in lei prevents any long range 
ordering, while stable ordered state is not forbidden in three-dimensional systems. 
In the intermediate case, due to the logarithmic diverging behaviour of the intergral 
it is less obvious to conclude and more refined analysis is required. In his famous 
book on phase transitions, Cardy uses a simplified version of Peierls argument on 
the existence of a phase transition at finite temperature in 2d in the case of discrete 
symmetry and extends the argument to continuous symmetry, showing that ordered 
ground state is unstable with respect to thermal fluctuations in this latter situation. 
Consider a spin system with only nearest-neighbour interactions — JSjSj and assume 
that the spins are represented by classical n— component vectors. An ordered ground 
state may be stabilized e.g. by symmetry breaking fields at some boundaries of the 
system. The variation of internal energy when a droplet of typical size / with spins 
progressively tilted in such a way that at the center of the droplet the spins are 
pointing opposite to the direction of the field is of order of 0(l 2 ) x J\S\ 2 (n/l) 2 where 
TV /I is the nearest neighbour spin disorientation. This result follows from integration 
over the droplet volume 0(l 2 ). Considering that the entropy is measured by the 
number of possible closed loops of size 0(1) in 2d, the entropy of the droplet is 
estimated as fc^ln//' where \i < z — 1 (z is the coordination of the lattice), so 
that eventually the free energy variation is AF(l) = n 2 J\S\ 2 — k B Tl\nfi. At any 
non zero temperature, increasing the size / of the droplet stabilizes the system and 
a spontaneously ordered ground state is thus impossible. In the case of discrete 



2 



Nematic phase transitions in two-dimensional systems 



symmetry, e.g. Ising model, the energy balance would be associated to the interface 
only and we would get AF(l) = 2J|S| 2 / — k bT I ha. fj,, showing that a transition to a 
stable ordered ground state is expected at a temperature around ksT c = 0(2J|S| 2 ). 
The two examples are illustrated in Fig. 1. In the latter situation, a conventional 




Figure 1. Evaluation of the free energy of a disordering droplet in two dimensions 
in the case of continuous symmetry (left) and discrete symmetry (right). 

phase transition toward an ordered phase is expected at finite temperature, while 
such an ordered phase may only be encountered at zero temperature in the first 
example [1,2]. On the other hand, an unconventional phase transition toward a 
quasi-long-range ordered state may take place at finite temperature, as we discuss 
in the next section. 



2. Two-dimensional electrodynamics and the XY model 

The celebrated XY model, 

H = ~JJ2 cos ( r ~ r>) ^ const + & ~ ^ r ') 2 ' ( 3 ) 

(r,r') ^ (r,r') 

admits a phenomenological description in terms of Coulomb gas. This is a well- 
known description, first given by Berezinskh and Kosterlitz and Thouless [3-5] and 
it is worth reminding its essential steps here. Before considering this model, imagine 
a point-charge q located at site r in a two-dimensional space, p(r) = qS(r — r ). The 
corresponding Coulomb potential, solution of Poisson equation V 2 0(r) = — Eq V( r ) 
might be written 

0(r-r o ) = --^ln^-^ (4) 

where a is chosen such that 0(a) = 0. Consider now a spin system living in two- 
dimensional space and let call u = V6 1 the distorsion field, where #(r) is the phase 
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field defined by S(r) = (cos 0(r), sin 0(r)). Due to the periodicity of the phase field, 
the distorsion field should obey the following relation § C (r Q ) u dl = 2n x integer 
where the contour integral is taken along a counterclockwise closed path around the 
point r . Using Stokes theorem, we may also write 

* udl — / d 2 rz(V x u) = 2ty x integer (5) 

/c(r ) J 

which implies 

V x u = 27r x integer x S(r - r )z. (6) 

The general solution of this equation consists of two terms, u = V^> — V x (z0). 
The curl of the first term being identically zero, we are led to a Poisson equation for 
the singular term V 2 = 2n x integer x S(r — rn) where the integer = n plays the 
role of the total (topological) charge enclosed by the contour C. From an energetic 
point of view, we may consider the kinetic energy 

PH= l -K j d 2 r |u| 2 (7) 

which, up to an inessential constant, is the continuous approximation of equation (3) 
with K = f3. J. After decomposition in ip and <fi parts, the cross-term vanishes and 
we get two independent contributions, 

[3H=±Kjd 2 r\V^\ 2 + ^jd 2 r|V0| 2 . (g) 
spin waves defects 

This is a fundamental relation which is the basis of the factorization of the par- 
tition function in a spin wave contribution and a defect (vortex) contribution, 
Z = Trexp(— j3H) = ZswZy. In the lattice version, the spin wave contribution 
to the Hamiltonian PHsw = |-K'Z)(r,r')(V'r — Vv) 2 is quadratic in Fourier space, 

^ sw = ^EI < ilVlV'q| 2 (9) 

A q 

(/i is the lattice spacing and ^ q the Fourier component of ip r ). This expression leads 
to the Gaussian model which implies [6] that 

(S ri • S r2 ) = e -i«^i-^) 2 > ~ |n - v 2 \- 1 ' 7fKK (io) 

where a temperature-dependent spin-spin critical exponent is found, r/sw — 
The low temperature (LT) phase of the XY model is a quasi-long-range ordered 
phase (QLRO), or a critical phase. When the temperature increases, vortices bounded 
in pairs produce further disordering of the system and a faster than linear increase of 
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the exponent rj with temperature up to the temperature where the transition takes 
place. The role of vortices is understood perturbatively through the calculation of the 
effective screened interaction energy between the topological charges. We consider 
a neutral Coulomb gas. Omitting a diverging contribution 0{\n.R/a) x J2i n i which 
would occur otherwise (i.e. for a non- neutral Coulomb gas), the vortices contribution 
to the Hamiltonian reads as 



I3H V = ~K J d 2 r<j)V 2 < 1 



r, ; - r 



= -2ttK ^Tn^ In ^ 3j - = '£ l PV(r i -r j ) (11) 



where V{r i — rj) is the Coulomb interaction between the charges. In the perturbative 
approach, we consider the effective interaction between charges rii and rij = —rii in 
the presence of another screening dipole. The following result comes out [7] 

e - / 3V eff (r i -r,) = e -^ (ri - rj )^ + congt xy 2 x J drr \-PV{r) + (){y*)), (12) 

where yo is the fugacity of a charge. The correction term y$ J dr r 3 ~ 2nK diverges at 
the deconfining transition 2nK c = 4 where the pairs of charges break. The presence 
of vortices increases the disordering of the system which, below T c , flows under 
renormalization to the zero-fugacity limit where the spin-wave limit is recovered 
with a renormalized temperature [5]. At the transition, the correlation function 
exponent takes a universal value 

"< = 2k = r < 13) 

The Heisenberg model (0(3)) on the other hand has no transition at any finite 
temperature (asymptotic freedom). An intuitive argument called escape to the third 
dimension is often reported. This is the observation than vortices cannot be stable 
for n > 3, since an infinitesimal amount of energy is able to produce a spin- wave- 
type excitation which eliminates the localized defect. Also a renormalization scheme 
was proposed by Polyakov to treat the non linear a— model close to the lower critical 
dimension d = 2 [8-11]. The flow of the coupling constant g (the temperature) under 
a change of scale s is governed by the (3— function, (3(g, d) = = (d — 2)g — (n — 
2)g 2 + 0(g 3 ) which shows that the model is disordered at any temperature when 
n > 2 in two dimensions, since the coupling always decreases and flows to zero under 
renormalization. 
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3. Nematics 

3.1. Definition of the model and of the observables 

We now come back to liquid crystals, the molecules of which may be idealized as 
long neutral rigid rods. They are likely to interact through electrostatic interactions, 
therefore Legendre polynomials appear for the description of the orientational transi- 
tion between a disordered isotropic high temperature phase and an ordered nematic 
phase. Lattice models of nematic-isotropic transitions capture the essentials of this 
extremely simplified description. The molecules are represented by n— component 
unit vectors S r (also called "spins"), located here on the sites r of a square lat- 
tice. The interaction between molecules is restricted to the nearest neighbour pairs 
(r, r'), the radial dependence being kept constant and the angular dependence enter- 
ing through a k— th order Legendre polynomial 1 , Pfc(S r • S r /), in terms of the scalar 
product between S r and S r /, S r -S r /. A coupling parameter J measures the interaction 
intensity. One obtains the following Hamiltonian of a lattice liquid crystal, 

ff = -iEft(Sr'Sr'). (14) 
(r,r') 

When the value of k in Eq. (14) is varied, new features may be expected, as in the 
case of symmetry-breaking magnetic fields hk cos k9 added to the XY model which 
change the phase diagram as investigated by Jose et al [12,7]. When the polynomial 
order k increases, one may expect a qualitative change in the nature of the transition, 
like in the case of discrete spin symmetries (Potts model) [13,14]. The value k = 2 was 
intensively studied. It still corresponds to the XY model for 0(2) spin symmetry, 
while it leads to the RP 2 or Lebwohl-Lasher model [15] for 3— component spin 
vectors. The nature of the transition in this latter case is still under discussion: an 
early study of Kunz and Zumbach [16] reported numerical evidences in favour of 
a topological transition, but more recently, several authors argued in favour of the 
absence of any finite-temperature phase transition [17,18], like in the Heisenberg 
case. Our own previous contributions [19,20] support the first scenario with QLRO 
at low temperature like in the XY model, but one cannot exclude a finite - but 
extremely large - correlation length which exceeds the maximum size available in 
numerical simulations. A recent study reported extremely convincing new evidences 
in favour of a topological transition [21], the transition being driven by topologically 
stable point defects known as | disclination points. Eventually, in the large— n limit, 
there is a proof of asymptotic freedom for values of k (in the interaction term 
(1 + cos9) k ) which do not exceed a critical k c ~ 4.537... [18]. Above this value the 
transition becomes of first order, a result which does not violate Mermin-Wagner- 
Hohenberg theorem, since the correlation length is finite at the transition. For finite 

1 Even order Legendre polynomials guarantee the local Z 2 symmetry S r — > — S r . 



6 



Nematic phase transitions in two-dimensional systems 



value of n, the question of the nature of the transition at high k is still a challenging 
problem, although there is a rigourous proof that the transition becomes of first-order 
for large enough values of k for arbitrary n > 2 [22,23]. This observation suggests 
to inspect the effect of a higher value of k also for the 0(2) model. In the context of 
orientational transitions in liquid crystals, Legendre polynomials rather than cos fc 
interactions are introduced, and we are led to the Hamiltonian of Eq. (14). 

In this report, we consider the behaviour of an Abelian spin model (0(2) rota- 
tion group) with P 4 — like spin interactions. We will refer to it as the P 4 0(2) model 
for simplicity. For 2— component vectors in a disordered phase, (cos 2 6) — | and 
(cos 4 #) = |. In order to keep the same symmetry in the interaction than in the 
P4 0(3) model already considered in the literature [24,25], but to normalise it be- 
tween and 1 in the limits of completely disordered and completely ordered phases 
respectively, we modify slightly the Hamiltonian, considering pair interactions of the 
formO^x) = AP 4 (a;)-|-const = ^(35x 4 — 30x 2 + ^). The corresponding Hamiltonian 
is thus defined by 

Hp 4 0(2) = —J ^2 Qi(Sr ■ S r /), (15) 

(r,r') 

with S r = (S%, S%), |S r | = 1. A qualitative description of the transition is provided 
by the temperature behaviour of the energy density, the specific heat, the order 
parameter and the susceptibility. The internal energy is defined from the thermal 
average of the Hamiltonian density, u(T) = (dL d )~ 1 (H) and the specific heat follows 
from fluctuation dissipation theorem, C V (T) = (L d T 2 ) _1 ((iJ 2 ) — (H) 2 ). Brackets 
denote the thermal average. The definition of the scalar order parameter (sometimes 
called nematisation) is deduced from the local second-rank order parameter tensor, 
M a/3 (r) = S?S% - \5 a P. After space averaging, the traceless tensor L' d £ r M Q/3 (r) 
admits two opposite eigenvalues ±|?7 corresponding to eigenvectors n + and n_. 
The order parameter density is defined after thermal averaging by M 2 (T) = (77). 
The associated susceptibility is defined by the fluctuations of the order parameter 
density, XM 2 (T) = 0((^ 2 )-(r ? ) 2 ). 

Simulations are performed using a standard Wolff algorithm adapted to the ex- 
pression of the nearest neighbour interaction [25,26]. The spins are located on the 
vertices of a simple square lattice of size L 2 with periodic boundary conditions in the 
two directions. Usually 10 6 equilibrium steps were used (measured as the number of 
flipped Wolff clusters) and 10 6 Monte Carlo steps for the evaluation of thermal aver- 
ages (the autocorrelation time at k B T / J = 0.2, L = 16 is of order of 30 MCS, hence 
the numbers of MC steps corresponds roughly to 3.10 4 independent measurements 
for the smallest size). A first qualitative description of the behaviour of the system 
is provided by the temperature dependence of thermodynamic quantities [27]. The 
specific heat has a maximum which does not seem to increase substantially with 
the system size. This might be the sign of an essential singularity around a temper- 
ature kpTc/ J ~ 0.70. From the order parameter variation, a smooth transition is 
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Figure 2. Energy, specific heat, order parameter M2 and corresponding suscepti- 
bility xm 2 vs T for the P4 0(2) model. The full lines are only guides for the eyes. 
The values of ks and J have been fixed to unity. 



suspected, since there is no evolution toward a sharp jump. The susceptibility dis- 
plays a non conventional behaviour at low temperature, increasing with the system 
size, which indicates a possible topological transition with a critical low temperature 
phase where the susceptibility diverges at any temperature. 

3.2. Characterization of the low-temperature phase 

We assume the existence of a critical phase at low temperatures as suggested by 
the temperature dependence results. The properties of the phase transition may be 
studied using 

i) Standard Finite-Size Scaling technique (FSS): in the critical low temperature 
phase of a model which displays a topological transition, the physical quantities 
behave like at criticality for a second-order phase transition, with power law 
behaviours of the system size. The difference is that in the critical phase, the 
critical exponents depend on the temperature and for any temperature below 
the transition one has e.g. 

M 2 (T) ~ L~^m 2 (t) ( 16 ) 
Xm 2 (T) ~ L^m 2 {t\ 
Here r] M2 (T) denotes the correlation function critical exponent, 

(Q 2 (cos(# ri - 0^))) ~ | ri - r 2 |-^ T ). (18) 
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ii) Rescaling of the density profiles (or "Finite Shape Scaling", FShS): conformally 
covariant density profiles or correlation functions are exepected at any tem- 
perature below the transition T c . They transform according to G(wi,w 2 ) = 
\w'(zi)\~ Xcr \w'(z 2 )\~ Xc 'G(zi, z 2 ) through conformal mapping w(z) where w la- 
bels the lattice sites in the transformed geometry (the one where the com- 
putations are really performed), while z refers to the infinite plane where 
the two-point correlations take the standard power-law expression G(z±, z 2 ) ~ 
I z\ — Z2\~ r,cr ), and x a = \r] a is the scaling dimension associated to the scaling 
field under consideration. Rather than two-point correlation functions, it is 
even more convenient to work with density profiles m(w) in a finite system 
with symmetry breaking fields along some surfaces in order to induce a non- 
vanishing local order parameter in the bulk [28-30]. The density m(w) will be 
M 2 (r) = (Q 2 (S r • h 9A )). In the case of a square lattice A of size L x L, with 
fixed boundary conditions along the four edges <9A, one expects (details may 
be found e.g. in [29]) 



m(w) ~ [ K (^)] 



k(w) = 3m 



\ 2Kw- 






sn — - — 


X 




Lj 







(19) 



Another conformal mapping which has been applied to many two-dimensional 
critical systems is the logarithmic transformation w(z) = In z = In p+i^ 
which maps the infinite plane onto an infinitely long cylinder of perimeter 
L. The correlation functions along the axis of the cylinder (let say in terms 
of the variable u = j-hvp) decay exponentially at criticality, G(ui,U2) ~ 
exp[—(u 2 —Ui)/£\, the correlation length amplitude on the strip being universal, 
£ = ^- [31]. Simulations at different temperatures in a system of size 10 x 10000 
were performed and the correlation function exponent thus follows from the 
linear behaviour 

ln(Q 2 [cos(0 U2 - e ui \) = const - ^^(u 2 - u ± ). (20) 

The resulting 77— exponent deduced from these different techniques is shown in 
Fig. 3. The variation of rj with the temperature is very similar to the one observed 
in the case of the XY model and confirms the QLRO nature of the LT phase of 
the model. The low temperature linear variation of this exponent is easily under- 
stood within the spin wave approximation. For 0(2) model with nearest neighbour 
interactions described by arbitrary polynomial in S r • S r /, one is led to an effective 
harmonic Hamiltonian ^JJ2(r,r') — @r') 2 - It yields power-law correlations, 

(cosm(# ri - e T2 )) = e-^-i-^) 2 ) ~ |n - r 2 p™< (21) 
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Figure 3. Temperature variation of the correlation function exponent rjM 2 (T) de- 
duced from conformal rescaling (open symbols) and FSS (filled symbols). Open 
triangles and diamonds, which correspond the largest systems seem quite reli- 
able. The dashed line shows the result of the spin- wave approximation t}m 2 (T) = 



11 k B T 
64 tt.J 



with a decay exponent given by 

2 2 

m m 

VrrU = -VXY = ^ (22) 

The comparison is made visible in the figure (with m = 2 and I = 128/11), and as 
expected, the lower the temperature the better the SW approximation. 



3.3. Critical behaviour at the deconfining transition 

Not only the low temperature behaviour of rj is interesting, the precise value of 
the f] exponent at the BKT transition where some deconfining mechanism should 
lead to the proliferation of unbinded topological defects is also of interest, since it 
really describes the universality class of the transition. An accurate value of the 
transition temperature is first needed. We performed a study of the crossing point 
of £7 4 Binder cumulant for very large statistics (30 x 10 6 MCS) and large system 
sizes (squares of L — 64, 80, 96 and 128 with periodic boundary conditions). The 
results shown in Fig. 4 indicate a transition temperature of /cbTbkt/</ = 0.7226. 

Then this temperature is used to perform Finite-Shape Scaling using the alge- 
braic decay of density profiles inside a square with fixed bounday conditions. These 
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Figure 4. Left: Crossing of the Binder parameter at the deconfining transition 
at a temperature IzbTbkt/J = 0.7226. Right: Size dependence of the correla- 
tion function exponent at Tbkt (L = 16, 32, 64, 128) and extrapolation to the 
thermodynamic limit. 



simulations are time-consuming, since the autocorrelation time increases in the low 
temperature phase when T evolves towards the deconfining transition and a rather 
large number of Monte Carlo steps is needed to get a satisfying number of indepen- 
dent measurements. For sizes L = 16, 32, 64, we used 10 6 MCS for thermalization 
and 30 • 10 6 for measurements, while "only" 20 • 10 6 for the largest size 128. The ex- 
ponential decay of two-point correlation functions along the torus cannot be applied 
at the BKT transition, since the system size being quite larger than in a square 
geometry, the number of MC iterations required is by far too large. In Fig. 4 we plot 
the "effective" exponent r) e s(L) measured at Tbkt for different system sizes as a 
function of the inverse size. An estimate of the thermodynamic limit value (L — ► oo) 
can be made using a polynomial fit (the results of quadratic and cubic fits are re- 
spectively 0.118 and 0.122), but it is safer to keep the three largest sizes available, 
L = 32, 64, and 128, for which a linear dependence of r] e g(L) with L~ l is observed. 
Taking into account the error bars, crossing the extreme straight lines leads to the 
following value for the correlation function exponent at the deconfining transition 

WTbkt) =0.122 ±0.007. (23) 
This value is essentially half the Kosterlitz value for the XY model. 

4. Summary and open questions 

The results obtained are essentially the following: 

- The P4 0(2) model displays a BKT-like transition with QLRO in the LT phase 
where SWA nicely fits the nematisation temperature-dependent exponents 
f](T) when T — > 0. 
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- The value of the exponent at the deconfining transition is supposed to reach a 
universal value. The numerical estimate is close to 0.125. 



We believe that our conclusions are safe in which concerns the existence of a 
QLRO phase at low temperature. The results may be understood through a naive 
comparison with clock model in 2d. Increasing the order of the interaction polyno- 
mial indeed increases the number of deep wells which stabilise the relative orienta- 
tion of neighbouring spins, and one is thus led to a system which is quite similar 
to a planar clock model with a finite number of states, unless the fact that here 
we keep a continuous spin symmetry which prevents from any "magnetic" long- 
range order at finite temperature. The clock model is in the Potts universality class 
when q = 3, but at q > 4, it displays a QLRO phase before conventional order- 
ing at lower temperatures [12]. Combining this analogy with the requirements of 
Mermin-Wagner-Hohenberg theorem for continuous spin symmetry gives a natural 
explanation for our results. For any type of nearest-neighbour interaction (in P 1 , 
P 2 or P 4 ) the behaviour seems to be always described by a BKT-like transition. 
The similar observation that a two-component nematic model renormalises in two 
dimensions towards the XY model was already reported in Ref. [32]. The transition 
is likely driven by a mechanism of condensation of defects, like in the XY model, 
but due to the local Z 2 symmetry not only usual vortices carrying a charge 1 are 
stable, but also disclination points carrying charges 1/2 should be stable. The role 
of these defects might be studied in a similar way than in the recent work of Dutta 
and Roy [21], by the comparison of the the transition in the pure model and in 
a modified version where a chemical potential is artificially introduced in order to 
control the presence of defects. 

Now, the observed value of the exponent r\ at the transition temperature seems 
a bit strange. A naive application of the mechanism discussed in section 2 leads 
to a value two times larger. Apart from minor modifications, the same procedure 
applies. The Hamiltonian (15) is of the form (7) with K replaced by Kl which does 
not affect the universal properties, but changes the critical temperature by the same 
factor. The other modification comes from the local Z 2 symmetry in the nematic 
model. This changes the charges n from integers to half-integers, which modifies the 
factor of 2ir in equation (11) to n/2. The divergence of the perturbative term in (12) 
then occurs at ^irKJ = 4, where the spin wave exponent (22) becomes i] c = m 2 /16 
which takes the value | when m — 2. We thus have a huge discrepancy of 100 % 
between the numerical determination and this prediction! A probable source of this 
discrepancy may be in the mechanism involved. There are no negative charge in the 
nematic model, so the pair interaction has a different structure and a more refined 
analysis is therefore desirable. 
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Note added in proofs 

We would like to thank S. Korshunov who draw our attention on a possible 
scenario explaining the discrepancy between the numerical value of rj and the pre- 
diction which follows from Kosterlitz-like arguments. Twenty years ago, he studied 
a similar planar XY model with mixed interactions [33]. In the parameters space 
of the problem, there is a zone where, starting from the high temperature phase, 
a BKT transition is first observed and then at lower temperature occurs a phase 
transition governed by the presence of solitons. This latter transition belongs to the 
two-dimensional Ising model (IM) universality class. According to this scenario, ours 
results would thus correspond to this IM transition, but the quantity that we study 
does not correspond to the order parameter of this transition (hence the value of 
rj is not the usual 1/4 of the IM universality class). QLRO would persist at higher 
temperatures, up to a limiting value where t]m 2 would indeed take its Kosterlitz-like 
value m 2 /16. 
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